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The morphological scaling properties of linear polymer films grown by vapor deposition polymer- 
ization (VDP) are studied by 1+1D Monte Carlo simulations. The model implements the basic 
processes of random angle ballistic deposition (F), free-monomer diffusion (D) and monomer ad- 
sorption along with the dynamical processes of polymer chain initiation, extension, and merger. The 
ratio G = D/F is found to have a strong influence on the polymer film morphology. Spatial and 
temporal behavior of kinetic roughening has been extensively studied using finite-length scaling and 
height-height correlations H(r,t). The scaling analysis has been performed within the no-overhang 
approximation and the scaling behaviors at local and global length scales were found to be very 
different. The global and local scaling exponents for morphological evolution have been evaluated 
for varying free-monomer diffusion by growing the films at G — 10, 10 2 , 10 3 , and 10 4 and fixing the 
deposition flux F. With an increase in G from 10 to 10 4 , the average growth exponent j3 « 0.50 was 
found to be invariant, whereas the global roughness exponent ot g decreased from 0.87(1) to 0.73(1) 
along with a corresponding decrease in the global dynamic exponent z g from 1.71(1) to 1.38(2). 
The global scaling exponents were observed to follow the dynamic scaling hypothesis, z g = a g /j3. 
With a similar increase in G however, the average local roughness exponent <xi remained close to 
0.46 and the anomalous growth exponent /9* decreased from 0.23(4) to 0.18(8). The interfaces dis- 
play anomalous scaling and multiscaling in the relevant height-height correlations. The variation 
of H(r,t) with deposition time t indicates non-stationary growth. A comparison has been made 
between the simulational findings and the experiments wherever applicable. 

PACS numbers: 82.20.Wt, 81.15.-z, 68.55.-a, 81.15.Gh 



I. INTRODUCTION 

Our motivation for gaining theoretical understanding 
of polymer thin film growth stems from their technolog- 
ical applications in microelectronic interconnects [l], , 
organic electronics [2|, and biotechnology. Various ex- 
perimental methods like vapor deposition polymerization 
(VDP) [l], UHl], ionization assisted polymer deposition 
p , sputtering growth Q , pulsed laser deposition || > 
and organic molecular beam deposition [10( have been 
developed to produce a variety of polymer thin films. 
Polymer film growth is complex compared to the con- 
ventional inorganic thin film growth process due to poly- 
mer's complicated structure and interactions that include 
internal degrees of freedom, limited bonding sites, chain- 
chain interactions, etc. Many experimental efforts have 
focussed on the formation of polymer thin films using 
VDP fllT - fl3l ]. In a typical VDP experiment, a wafer 
(2-D substrate) is exposed to one or more volatile gas 
phase precursors that produce free-monomers. The free- 
monomers impinge on the substrate at random locations 
and react on the substrate surface to produce the de- 
sired deposit. Polymer thin films grown by VDP are 



made up of long polymer chains formed through the poly- 
merization reaction occurring during the growth process. 
The polymerization process involves the interaction of 
two free-monomer molecules in a chemical reaction to 
initiate a dimer (polymer chain of length = 2). The free- 
monomers moving towards the substrate are consumed 
by either of the two processes: first being chain initia- 
tion in which new polymer molecules are generated; and 
secondly, chain propagation in which the existing poly- 
mer molecules are extended to higher molecular weight. 
Besides these two mechanisms, the free-monomer adsorp- 
tion, diffusion, and polymer merger can be considered as 
other mechanisms that determine the overall film mor- 
phology. The chemical nature of the linear polymer chain 
restricts the number of bonding sites. A free-monomer 
can only bond to either of the two active ends of a poly- 
mer, or to another free-monomer. This bonding con- 
straint leads to the formation of an entangled or an 
overhang configuration, which blocks the region it cov- 
ers from the access of other incoming free-monomers. In 
conventional physical vapor deposition (PVD) processes 



14| . atoms can nucleate at the nearest neighbors of the 
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nucleated sites and atomistic processes such as surface 
diffusion, edge diffusion, step barrier effect, etc. effect 
the growth, resulting in the films being compact and 
dense Recent investigations by Zhao et al. pjj 

have shown that the submonolayer growth behavior of 
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VDP is very different from that of PVD due to long 
chain confinement and limited bonding sites, indicating 
that the detailed molecular configuration can drastically 
change the growth behavior [2(|. In experiments, the 
growth behavior of polymer thin films have been investi- 
gated through their morphological evolution study. The 
VDP processes for producing Parylene-N (PA-N) films 
typically are far from equilibrium. The precursor ma- 
terial di-p-xylylene (dimer) is sublimed at 150°C and 
then pyrolized into free-monomers at 650°C. The free- 
monomers impinge at random angles onto the Si-wafer 
at room temperature and eventually condense and poly- 
merize to form the polymer film. By varying the growth 
rates in the PA-N growth experiments, Zhao et al. [2l[ 
reported an average roughness exponent a = 0.72 ± 0.05 
and an average growth exponent f3 = 0.25 ± 0.03. How- 
ever, by considering the tip effect of the atomic force 
microscope [22l |. the range of a was estimated to be be- 
tween 0.5 to 0.7 and the authors found the absence of 
dynamic scaling hypothesis in the PA-N film growth. In 
the recent experiments done by Lee et al. [231 ] . the au- 
thors observed unusual changes in the roughening be- 
havior during the poly(chloro-p-xylylene) growth. In the 
early rapid growth regime, they observed (3 = 0.65 (larger 
than the random deposition /? = 0.5) and upon com- 
plete coverage of the substrate (around d — 10nm), they 
found /? = 0.0 and the interface width did not evolve with 
the film thickness. Finally, during the continuous growth 
regime, the surface roughness again was found to increase 
steadily with a new power law of 8 = 0.18 [HI, which is 
close to the results of Zhao et al. [21j • Of the known the- 
oretical results for dynamic roughening, the MBE non- 
linear surface diffusion dynamics proposed by Lai and 
Das Sarma [24[ predicted similar exponents as obtained 
in the experimental study of Zhao, et al. [21]. However, 
their nonlinear surface diffusion theory could not explain 
the findings of varying local slope in the experiments. 
Zhao et al. proposed a stochastic growth model based on 
bulk diffusion 19], which correctly predicted the kinetic 
roughening phenomena observed in their experiment (2lj 
but lacks the details on how polymers evolve. Insufficient 
theoretical studies coupled with inconsistencies in the ex- 
periments motivate us to model the polymer film growth 
and seek a better understanding of the growth processes 
that determine the roughening mechanism. In this pa- 
per, we study the 1+1D lattice model for the polymer 
films grown by VDP and examine the effects of random 
angle deposition, free-monomer diffusion, free-monomer 
adsorption in determining the evolution of the film's mor- 
phology. 



II. MODEL AND METHOD 

In our simulation, the free-monomers were deposited 
at random angles on a 1-D substrate of lattice size L at 
a deposition rate F (in units of monomers per site per 
unit time) . The KISS random number generator (25| was 
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FIG. 1: Schematic of the 1+1D growth model. 1(a): Free- 
monomer deposition at random angles; 2(a): Adsorbed free- 
monomer diffuses along along the substrate; 2(b): Adsorbed 
free-monomer diffuses on a polymer chain; 3(a,b): Poly- 
mer chain initiation resulting from random angle deposition; 
3(c,d): Polymer chain initiation resulting from free-monomer 
diffusion on substrate and polymer chain respectively; 4(a,b): 
Free-monomer deposits onto the active end of a polymer chain 
resulting in chain propagation; 4(c) Chain propagation due to 
free-monomer diffusion. 



employed and one deposition time unit corresponded to 
the deposition of L free-monomers. The incident free- 
monomers were released from a height three lattice units 
above the highest point on the surface with an initial 
abscissa randomly chosen from 1 < x < L. The de- 
positing free-monomers have a uniform "launch angle" 
distribution which corresponds to a nonuniform flux dis- 
tribution J (6) ~ l/cos(t9) of particles above the surface, 
where is defined as the angle between the direction of 
impinging monomer and the substrate normal [26l [27] ]. 
Our VDP growth model is similar to the square-lattice 
disk model studied by Ref. (2(| along with additional 
constraints of free-monomer diffusion, limited bonding, 
polymer-initiation, propagation and merger. The par- 
ticles followed a ballistic trajectory until contacting the 
surface. The impinging free-monomer was then moved 
to the lattice position nearest to the point of contact. 
The deposited free-monomers were allowed to diffuse via 
nearest-neighbor hops with a diffusion rate D (nearest 
neighbor hops per monomer per unit time). A free- 
monomer that is deposited on top of an existing poly- 
mer chain gets adsorbed on the chain and is also allowed 
to diffuse. The excluded volume constraint was imple- 
mented by rejecting the diffusion or deposition moves to 
an already occupied site. As the monomer coverage in- 
creases on the substrate, the polymer film grows along 
a direction perpendicular to the substrate. This "two 
dimensional" growth is often referred to as the 1 + 1D 
growth in literature. Figure Q] shows the schematic of 
various processes that occur during the non-equilibrium 
film growth on a 1-D substrate of length L with peri- 
odic boundary conditions. Processes 1(a), 3(a), 3(b), 
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(a) G = 10, t = 15 



(d) G = 10, f = 15 



(b) G = 10, t = 30 





















(c) G = 10, t = 60 




(e) G = 10, f = 30 




FIG. 2: Snapshots of 
the polymer films grown 
on a L = 512 substrate. 
Figs. (a, b, c) were 
obtained for G = 10 
at deposition times t — 
15, 30, and 60 units re- 
spectively; whereas Figs, 
(d, e, f) were obtained 
for G = 10 6 and de- 
position times t — 15, 
30, and 60 units, respec- 
tively. One deposition 
time unit corresponds to 
a deposition of L free- 
monomers. For a clear 
view, the growth front 
(shown in red) is dis- 
placed vertically by 100 
pixels and the longest 
polymer chain is high- 
lighted in black. 



4(a), and 4(b) show the gas phase free-monomers de- 
positing onto the substrate at random locations with 
uniform launch angle distribution. These free-monomers 
get adsorbed either on the substrate (shown by process 
2(a)) or on the polymer chains (process 2(b)). Adsorbed 
free-monomers are allowed to diffuse along the adsorbent 
to any of the nearest-neighboring unoccupied sites with 
equal probability, the rate of diffusion D is assumed to be 
the same on both the substrate and the polymer chain. 
We define the chain length s of a polymer as the number 
of monomers forming a polymer chain. When an im- 
pinging free-monomer encounters another free-monomer 
on the substrate as its nearest neighbor (process 3(a)), 
both are frozen and undergo a chemical reaction to form 
a dimer (polymer chain of length s = 2); polymers with 
chain length s = 3 can also be formed after deposition 
(process 3(b)). Free-monomer diffusion along the sub- 
strate (process 3(c)) and along the polymer chain (pro- 
cess 3(d)) can also result in the polymer chain initiation 
128]. When an impinging free-monomer encounters the 
active end of a polymer chain, it attaches itself to the 
chain and increases the chain length by one unit (4(a), 
4(b)). A diffusing free-monomer can meet an active end 
of a polymer chain in its neighborhood and get bonded to 
that polymer chain (process 4(c)). In linear polymer sys- 
tems, the free-monomers are allowed to form a maximum 
of two chemical bonds and at any given time only the two 
ends of the polymer chain arc chemically active, result- 
ing in the chain propagation at these two end locations. 
The chain portion (of the polymer) excluding the two 
chemically active ends, is not allowed to form chemical 
bonds with neighboring free-monomers. Free-monomers 
can however be physically adsorbed on the chain and 



can diffuse along the chain (process 2(b)). Another in- 
teresting process that occurs during the film growth is 
the polymer-merger. Two different polymers can merge 
when their respective active-ends meet as nearest neigh- 
bors. During the polymer-merger, the nearest neighbor- 
ing ends of merging polymers react chemically to join the 
two polymers into one longer polymer chain with higher 
molecular weight. The resulting polymer chain is left 
with two active-ends, one from each of the parent poly- 
mers. In the case when both the active-ends belonging to 
the same chain appear as nearest neighbors, the chemical 
bond between the ends is prohibited. In this study we do 
not attempt to study the effects of re-emission which al- 
lows the free-monomers to "bounce around" before they 
settle at appropriate sites on the surface. Instead, we as- 
sume that the impinging free-monomers will always stick 
to the particle that comes on its deposition path (pro- 
cesses 4(a), 4(b)). At each stage of the simulation, 
either a deposition or a diffusion is performed. In order 
to keep track of the competing rates of diffusion and de- 
position we adapted the method suggested by Amar et 
al. [29j and carried out the deposition with a probability 

PF, 

PF = [i + aUg]' (1) 

where N\ is the free-monomer density (per site) and G = 
D/F. The diffusion was carried out with the probability 

PD, 

N 1 x G 

PD = 1 - pF= [l + N 1 xG\- (2) 
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In our simulations the incoming free-monomer flux F 
was fixed for different D, thus an increase in D was 
parametrized as an increase in G. Throughout the 
growth process the list of all free-monomers and poly- 
mer chains were continually updated. If a free-monomer 
encountered another free-monomer or an active end of a 
polymer as its nearest neighbor, it was added to the poly- 
mer chain and removed from the free-monomer list. In 
cases where a free-monomer was the nearest neighbor to 
the active ends of more than two polymers, we selected a 
random pair of polymers and performed polymer-merger. 



III. RESULTS 

A. Surface Morphology 

In Fig. [2] we show typical snapshots of the polymer 
films generated using L = 512 substrate for two extreme 
cases: G = 10 (Figs. Hi, b, c) and G = 10 6 (Figs. Hi, 
e, f) after a deposition time of t = 15 (Figs. [3h., d), t = 
30 (Figs. [2p, e), and t = 60 (Figs. [2fc, f) respectively. 
For both the values of G, the films show the presence of 
columnar structures, overhangs, and unoccupied regions. 
These structures were observed to persist throughout the 
growth process. Presence of these morphological struc- 
tures can be explained by the shadowing effect inherent in 
the growth process and is attributed to the cos^ distri- 
bution of the impinging free- monomer flux [13, [30l - l33| . 
Shadowing effects arise when the columnar structures 
of the surface "stick out" and shadow their neighbor- 
ing sites, thus inhibiting the growth in their neighboring 
sites. Due to the angular flux distribution of the imping- 
ing free-monomers the taller surface features prevent the 
incoming flux from entering the lower lying areas of the 
surface. 

For comparable deposition times t, the films grown 
at G = 10 (Figs. [5^, b, c) are characterized by small 
unoccupied-regions and short polymer chains, resulting 
in shorter, denser, and compact films. Whereas for G = 
10 6 (Figs. e, f), the films are characterized by large 
unoccupied-regions and longer polymer chains resulting 
in taller, more porous, and less dense films. For a rel- 
atively low diffusion rate at G = 10, the free-monomers 
deposited on the film have a higher probability of encoun- 
tering another impinging free-monomer as nearest neigh- 
bor and thereby initiating new polymers. Many such 
polymer initiations inhibit the occurrence of unoccupied- 
regions and make the film dense and compact. In con- 
trast, at a higher diffusion rate of G — 10 6 , the free- 
monomers have a higher probability to diffuse upward 
towards the growth front and the upward diffusion of free- 
monomers is favored due to the non-symmetric nature 
of the lattice potential associated with diffusion over a 
step (l5| . The diffusing free- monomers arrive towards the 
growth front and bind to the active ends of the polymers 
and increase their chain length. This explains the oc- 
currences of longer polymer chains at higher G observed 
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FIG. 3: Generalized correlation function H q (r,t) calculated 
for L — 2000 and t — 1000. The data are averaged over 
200 independent runs and the error-bars are smaller than the 
symbol size. 

in Fig. [2] (the longest polymer chains are highlighted 
in black). Throughout the growth process the longest 
chains for G = 10 6 are much longer than those obtained 
at G = 10. In Fig. [2] for a fixed t even though the films 
have the same number of particles, the film morphology 
looks different for G = 10 and G = 10 6 . The difference 
arises from the variation in the position of growth front 
h(x,t), defined as the set of occupied sites in the film 
that are highest in each column, and x represents the 
horizontal lattice site on the substrate. The growth front 
h(x,t) studied here is a crude approximation of a more 
complex aggregate that is growing. However, our method 
of quantifying the growth front is justified because, in the 
AFM experiments, the measured 1-D height-height cor- 
relation function is based only on height profiles along 
the fast scanning axis. Moreover the finite size of the 
AFM tip is known to distort the growth front and the 
measured growth front is a convolution in which the in- 
teraction with the tip dilates the surface details (22[. In 
Fig. [2] the growth fronts are shown vertically displaced 
by 100 pixels in the growth direction for clarity. The 
height fluctuation frequencies in h(x,t) are observed to 
occur at different length scales depending on the ratio 
G. For a fixed t, it is intuitive to think that the height 
profiles corresponding to G — 10 6 (Figs. [5U, e, f) arc 
more rougher than those obtained using G = 10 (Figs. 
[2^,, b, c). A detailed analysis of the dynamics of interface 
roughness evolution is presented in the later section. 

The gth-order generalized height-height correlation 
function H q (r, t) (defined below) is commonly used in 
identifying multi-affine surfaces [TH, HH . 

H q (r,t) = {(\h(x + r,t)-h(x,t) (3) 

The scaling properties of multi-affine surfaces can be de- 
scribed in terms of an infinite set of Hurst exponents 
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. 4: Variation of havg as a function of deposition time t 
for different G and L — 512. The data are averaged over 500 
independent runs. The statistical errors are smaller than the 
symbol size. 



FIG. 5: Variation of r(G) as a function of G. The data were 
obtained for L = 512 and averaged over 500 independent runs. 
The statistical errors are smaller than the symbol size. 



[IE Hl[ h q which are obtained using, 

H,(r,t)~rV (4) 

For multi-affine interfaces, the exponents h q are known 
to vary with q [l6j]. In Fig. [3l we plot the generalized 
correlation function H q {r, t) using q = 2, 3, 4, and 5 for G 
= 10 at deposition time t — 10 3 . For smaller r, we observe 
a power- law dependence of H q (r, t) on r in accordance to 
Eq. ((4]). The slopes of the log-log plots shown in Fig. [3] 
depend on q and indicate the presence of multi-scaling in 
the films generated by our VDP growth model. 

B. Average Height and Growth Rate 

The average height of the growth front h avg (t), is de- 
fined as, 

1 L 

h avg {t) = -jThfat), (5) 

X = l 

and quantifies the overall thickness of the film. Figure 
2] shows the plot of h avg versus t for different G and 
L = 512. For the studied values of G, the h avg versus t 
plots show a linear relationship. This linear behavior is 
a consequence of restricting the growth to 1-D and the 
excluded volume constraint implemented in the growth 
process. In Fig. [4] with an increase in G, we observe a 
systematic increase in the slope of the h avg versus t plot. 
We define the growth rate R of the polymer film as, 

R= dhav IL 

at 

In general one expects the polymer film's growth rate 
to be proportional to the incoming free-monomer flux 
only, i.e. R oc F and it is natural to expect R to be 
independent of G(= D/F) since F(= 1/L) is a constant 
for a fixed L. However, in Fig. 2] our simulations show a 



strong dependence of R on G as well, and the growth rate 
R is observed to increase monotonically with G. This G 
dependent growth rate R indicates that the growth rate 
is effected by monomer diffusion directly, i.e., there is a 
net uphill monomer diffusion current that is responsible 
for this R — G relationship. We thus incorporate an ad- 
ditional term r(G) that determines R in addition to its 
dependence on F. 

R(G) = R + r(G), (7) 

where R is the growth rate due to the random deposition 
flux only (D — 0) and r(G) is the growth rate induced 
only by the uphill diffusion of free-monomers. In Fig. [5] 
we show the variation of r(G) as a function of G. We 
find a monotonic increase in r(G) with an increase in G. 
This shows the strong influence of G in determining the 
polymer film's growth rate. Specifically, at higher G we 
find r(G) oc G 0162(3) . This indicates that the growth 
rate due to uphill diffusion r(G) asymptotically follows a 
power-law dependence with G. 

C. Characterization of the Kinetic Roughening 

The morphology of the growth front can be charac- 
terized by studying its spatial and temporal evolution. 
In morphological scaling studies, typically two kinds of 
scaling behaviors are associated with the roughenin g ki - 
netics: the global scaling and the local scaling [36m 8j |. 
The growth models with anomalous kinetic roughening 
are known to have different scaling exponents in the lo- 
cal and the global scales. The local roughness exponents 
have been employed in studying the irregularly growing 
mound mor phology and are often used in the experimen- 
tal analysis [2l|, |39( . In general, the local roughness expo- 
nent cti and the global roughness exponent a g take differ- 
ent values depending on the type of scaling exhibited by 
the growth process. In the case of super-rough surfaces 
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FIG. 6: Lateral film density profiles p (y) of polymer films 
grown on a L = 512 substrate with G = 10 at deposition 
times of t = 10, 20, 40, and 60 respectively. The statistical 
errors are smaller than the symbol size. 



generated by nonequilibrium MBE growth models [40| 
and growth models with horizontal diffusion [4l[ , the as- 
sumption of the equivalence between the global and local 
descriptions of the surface is not valid and such behav- 
ior has been termed as anomalous scaling [IE E3 . The 
differences in the global and local scaling exponents have 
been attributed to the super-roughening and intrinsically 
anomalous spectrum observed in the anomalous scaling 
of surfaces [43]. For studying the morphological evo- 
lution of VDP generated polymer films, it is essential to 
identify their steady growth regime. To do so we em- 
ployed the lateral film density p (y) (at a height y lattice 
units) of polymer film defined as, 



p(y) 



N(y) 



(8) 



where N(y) represents the number of occupied lattice 
sites at a height y above the substrate. Calculations were 
performed on films with L = 512, G = 10 at deposition 
times of t — 10, 20, 40, and 60, respectively. Figure [5] 
shows the variation of p(y) with t. Similar plots were 
obtained for other G values also. Three distinct regions 
of the polymer film growth referred to as the interface 
region, the steady-state film region and the growth front 
were identified based on the Ref. [2(j- To evaluate the 
scaling exponents for kinetic roughening, we only con- 
centrate on films with the deposition time t 3> ti nter face- 



1. Global Scaling Behavior 

The global scaling of the growth front can be deter- 
mined by studying the lattice size dependent interface 
width w(L,t) defined as, 



w(L,t) = 



\ 



1 L 

J^iHi) - h avgY 



(9) 



In many 1-D morphological growth processes such as bal- 
listic deposition [44], |45| , Eden model [45J , and solid-on- 
solid models [46j], w(L,t) usually follows the scaling law, 



w(L,t) ~ t 







(t < t x ), 



(10) 



where exponent ft is known as the growth exponent that 
characterizes the time-dependence of surface roughening. 
For any given L, the power-law increase in w{L, t) does 
not continue indefinitely with t, and is followed by a sat- 
uration regime during which the w(L,t) reaches a satu- 
ration value w sa t ■ The power-law growth regime and the 
saturation regime are separated by a crossover time tx- 
As L increases, w sa t also follows a power-law (l5j . 



w(L,t) 



L a <> 



(t > t x ), 



(11) 



and a g is referred to as the global roughness exponent 
[l5j . The time tx at which the behavior of w(L, t) crosses 
over from Eq. (fTU)) to that of Eq. (fTTj) depends on L and 
scales as, 



t x ~U- 



(12) 



where z g is called the global dynamic exponent. Typi- 
cally the scaling exponents are independent of specific 
interactions involved in the growth process and depend 
on the dimensionality and symmetries of the system 
[IE EE EE S3- F° r some growth processes, the expo- 
nents a g , 0, and z g are unified using the dynamic scaling 
hypothesis [45[ (also known as the Family- Vicsek scaling) 



given as, 



w(L,t) ~ L a ^(t/L z <>) 



(13) 



where ^ (t/ ' L z <>) is referred to as the scaling function and 
satisfies, 



const 



(x < 1) 
(>»1), 



and 



Otg/0. 



(14) 



(15) 



In Fig. [7] we show two representative plots of w(L, t) 
versus tona log-log scale for (a) G = 10, (b) G = 10 3 and 
varying L. From both Figs. [TJa) andJTJb), for t <€. tx, 
the log{w{L,t)) versus log(t) shows a linear dependence 
implying a power-law behavior. For the same G and 
t <C <X) the w(L,t)-t plots overlap with one another for 
different L. We estimated (3 by fitting Eq. (fTU|) to the 
plots of w(L,t) in Fig. [7] for the film region t tx- 
We obtained average (3 = 0.509(8) and 0.51(2) for G = 
10 and 10 3 respectively. For other G values, the (3 ob- 
tained was close to 0.50 and we observed an invariance of 
j3 (within the error-bars) with G. The statistical average 
and error bars in f3 were obtained from 1.8 x 10 3 inde- 
pendent simulations. The VDP model studied here is 
similar to the ballistic deposition model with additional 
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FIG. 7: Time variation of the interface width w(L, t) for (a) G = 10 and (b) G = 10 3 with L = 100, 200, and 300. The data 
are averaged over 1.8 x 10 3 independent runs. 
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FIG. 8: Variation of w sa t with L for G varying from 10 to 
10 4 . The data are averaged over 1.8 x 10 3 independent runs 
and the statistical errors are smaller than the symbol size. The 
error in the exponents were obtained from the curve- fit. 



degree of freedom including diffusion, polymer initiation, 
extension and merger. However, unlike the atomic diffu- 
sion in MBE growth, the free-monomer diffusion is con- 
fined by the linear geometry of the polymer chain. From 
Fig. [2] one can observe that in most cases the polymer 
chains are more or less perpendicular to the substrate. 
Since the free-monomers can only move along the poly- 
mer chain, most of the diffusion happens in vertical direc- 
tion rather than in lateral direction (which is the case for 
MBE growth). Yet, the vertical diffusion does not con- 
tribute significantly to extra roughness increasing or de- 
creasing (the total particle number should be conserved) 
in the growth front due to the porous nature of the film 
as compared to the random deposition model. Thus, it is 
expected that the growth exponent will be close to that 
of random deposition (/3 ~ 0.50). We calculated w sat 



by averaging w(L, t) for t 3> tx from the data shown 
in Figs. [7][a) andJTfb). Figure [5] shows the plot of w sat 
versus L on the log-log scale for G varying from 10 to 
10 4 . The exponent a g was estimated for each G using 
Eq. (fTTj) . With an increase in G from 10 to 10 4 , a g was 
found to decrease from 0.87(1) to 0.73(1), the error-bars 
in a g were obtained from the curve fitting. The exponent 
a g is known to be closely related to the surface fractal 
dimension [l5[ . The smaller the a g , the larger is the frac- 
tal dimension. Our observation of a decrease in a g with 
an increase in G shows that the fractal dimension of the 
growth front increases with G, i.e. there are more spatial 
frequency fluctuations in the film morphology with an 
increase in G. This finding is consistent with the growth 
front profiles shown in Fig. [5J This result also demon- 
strates that G induces a large effective vertical growth 
rate R(G) (shown in Fig. [SJ and the large R(G) in turn 
produces a much rougher film surface. 

To determine whether the VDP process obeys the dy- 
namic scaling behavior, we start with an assumption that 
VDP growth follows the dynamic scaling hypothesis and 
obtain z g through Eq. (fT"5"|) . As w(L,t) scales with both 
t (Eq. ((TD])) and L (Eq. dTTJ)) we can rescale the w(L,t) 
curve shown in Fig. [7] by plotting w/L ag versus t/L Zg 
to see whether those curves "collapse". For G = 10 we 
obtain a g = 0.87(1), /3 = 0.509(8) and according to Eq. 
(HU z g = 1.71(1) and for G = 10 3 we get a g = 0.780(8), 
P = 0.51(2), and z g = 1.53(2). We use the data from 
Figs. E(a) andEKb) and divide w(L,t) by L Q ». This 
shifts the curves of varying L vertically on the log-log 
scale. According to Eq. (fTTj) . these curves now satu- 
rate at the same value of the ordinate w/L ag , however 
their saturation times do not overlap. We then rescale 
the time axis and plot t/L Zs for both cases of G. This 
rescaling of time axis according to Eq. (|T2"j) leads to a 
horizontal shift of the curves and the curves now satu- 
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FIG. 10: Auto-correlation function C(r, t) calculated for L — 
2000 and t — 1000 (averaged over 200 independent runs). 

rate at the same abscissa t/L Zg . In Figs. HJa) andHJb) 
we show the rescaled plots of w(L, t) for G = 10, 10 3 and 
consequently observe the "collapse" of individual curves 
for varying L onto a single curve. This characteristic 
"collapsed" curve shown in Figs. M, a ) and Mjj) is the 
scaling function ty(t/L Zg ) mentioned in the Eq. (|14p. 
The scaling functions obtained in Figs. Uta) andEJb) are 
observed to follow Eq. (|14[) for both cases of G = 10 and 
10 3 . For other G values, we observed similar "collapse" 
behavior indicating the global dynamic scaling of growth 
fronts of the polymer films grown using VDP. 



2. Local Scaling Behavior 

The local scaling behavior of the growth front can be 
understood from studying the spatial correlation func- 
tions: the auto-correlation function C{r,t) and height- 
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FIG. 11: Estimates for column separation A and lateral 
correlation length £ obtained for G = 10(circles) and G = 
10 3 (triangles) . Open symbols correspond to £ and filled sym- 
bols correspond to A. The data are averaged over 200 inde- 
pendent simulations with L — 2000 and t = 1000. The error 
in the exponents were obtained from the curve-fit. 



height correlation function H(r, t) that describe the local 
properties of growing interfaces, 



C(r,t) = (h(x + r,t) ■ h(x,t)}, 



(16) 



H(r,t) = {([h(x + ri t)-h(x,t)} 2 )} = 2[w 2 -C(r,t)}, 

where r is the translation distance also referred to as the 
lag or slip [48j and (...) denotes a spatial average over 
the entire system. The functions C(r,t) and H(r,t) are 
directly related as shown in Eq. (JT7J) and differ only by a 
constant pre- factor of 2w 2 . The H(r, t) scales in the same 
way as the interface width and is often used in studying 
the kinetic roughening [l5| . In order to obtain accurate 
parameters from the correlation functions, it is impor- 
tant to account for the accuracy of statistical averages. 
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FIG. 12: Height-height correlation lunction H(r,t) for (a): G = 10, (b): G = 10 3 using L = 2000 at t = 500, 750, and 1000. 
The data are averaged over 200 independent runs and error in the exponents were obtained from the curve- fit. 



For the random Gaussian model surface studied by Ref. 
(48j, converging H(r, t) were obtained within an order of 
{£/L) d ^ 2 where d is the system dimension and £ is the lat- 
eral correlation length. The accuracy depends not only 
on the number of data points, but also on the sample size. 
It is the ratio £/L, not the number of data points, that 
determines the accuracy. Once the ratio £/L is known, 
one may not be able to increase the accuracy no mat- 
ter how many data points are collected [48| . This rule is 
different from the law of large numbers for independent 
random variables. This important difference needs to be 
recognized while studying spatially correlated systems. 
Ideally, one would like to have \J£JL <C 1, i.e, L — > oo. 
However, due to the computational constraints we per- 
formed calculations using L = 2000 and t — 1000. For 
a random self-affine surface, C{r, t) usually decays to zero 
with an increase in r. The shape of the decay depends 
on the type of the surface and the decay rate depends 
on the distance over which two points x, x + r become 
uncorrelated. In Fig. [2] the polymer growth front does 
not appear to be a random rough surface, instead it has 
regular fluctuations of the columnar structures. Figure 
[TU] shows the C(r, t) for L = 2000 after a deposition time 
t = 1000 for G = 10, 10 2 , and 10 3 used to characterize 
the surface morphology as a function of t. We can define 
two different lateral length scales: the lateral correlation 
length £ and the average column-separation A. The lat- 
eral correlation length £ defines a representative lateral 
dimension of a rough surface and is estimated through 
C(r,t) using C(r = £, t) = C(0,t)/e. Within a distance 
of £ the surface heights of any two points are correlated. 
The parameter A characterizes a wavelength-selection of 
a surface and is determined by measuring the value of 
r corresponding to the first zero-crossing of C(r,t) [48| . 
The variation of A with t represents how the columnar 
structures coarsen with deposition time. In general, the 
evolution of the columnar feature size follows a power law 



with t given by [26[ 

A oc t p , (18) 

and p can be referred to as the coarsening exponent. 

In Fig. QHwe plot A versus t for L = 2000 and G = 10 
and 10 3 along with the estimates for p. With an increase 
in G from 10 to 10 3 , p is found to remain close to 0.57. 
The invariance of p indicates that the coarsening of the 
columnar structures follow a power law that is unaffected 
by the ratio G. And we believe that this is dominated by 
the shadowing effect due to monomer vapor coming from 
all different angles uniformly. 

The local roughness exponent ai of the growth front 
can be obtained from H (r, t) using flU , [36j 

H(r,t)~r 2a ' (r«£)- ( 19 ) 

In Figs. H2Ta) andH^b) we plot H(r, t) for t = 500, 750, 
and 1000 along with the estimates for ai for L — 2000 
and G — 10, 10 3 respectively. For a given G and varying 
t, the estimates of a/ were observed to remain invariant 
within statistical-errors. From Figs. [TSTa) and [12T b). 
we obtained an average a/ of 0.470(3) and 0.460(3) for 
G = 10 and 10 3 respectively. A comparison between a g 
and ai (shown in Fig. [5] and Fig. IT2"j) shows a g > ai 
for studied G and the scaling behavior is observed to be 
different at short and large length scales. For small r, the 
plots of H(r, t) do not overlap at varying coverage and 
show the presence of non-stationary anomalous scaling 
p^ |. The vertical temporal shift in the H(r,t) observed 
in Fig [12] is due to the difference between a g and u\ and 
indicates the presence of anomalous scaling. The mech- 
anisms that lead to anomalous scaling can be separated 
into two classes: super- roug hening (a g > 1) and intrinsic 
anomalous scaling [38, 49', |50|. The anomalous growth 
exponent /?„ = (a g — cti)/z g [43j measures the difference 
between a g , ai and can be obtained from the scaling of 
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FIG. 13: Plots of rescaled H(r,t) showing the data collapse for (a): G = 10, and (b): G = 10 3 using L = 2000 at t = 500, 
750, and 1000. The data are averaged over 200 independent runs and error in the exponents were obtained from the curve-fit. 
The straight lines are plotted as a guide to the eye. The scaling agrees with Eq. (|21|l . 



H(r,t) [5C 



H(r,t) = r 2a °g A (r/t 1 ^), 



(20) 



and the anomalous scaling function <?a(u) (43T. |50I| satis- 
fies 



9a{u) 



(«« 1) 
(«»1) 



K\ = 2(a g - ai) 



k 2 = 2a„ 



(21) 



(22) 



(23) 



In Figs. [Ufa) and[T3Tb) we show the plots of H(r, t)/r 2a <> 
versus r/t 1 / 7 -" and obtain the "data-collapse" for G = 10 
and 10 3 respectively. The "collapsed" curves shown in 
Figs. ITBT a) and fTST b) are the scaling functions fiu(u) (for 
G = 10 and 10 3 ) mentioned in Eq. ([20]). For both G, the 
scaling functions obtained in Fig. [13] satisfy Eq. (12T1) in 
accordance with the theory for anomalous scaling [43, 50] . 
The exponents k% and k 2 were obtained from qa{u) plots 
using Eq. flST}. For G = 10 we have k x = 0.823(8), n 2 = 
1.731(2) and for G = 10 3 we obtained m = 0.621(4), 
«2 = 1.555(8). For both G, we find that the exponents K\ 
and K2 obtained from the curve-fit satisfy Eqs. (|22l) and 



for the numerical estimates of a g and a; obtained 



from Eq. (p} and Eq. ([T9]) . 

The dynamic scaling exponents of the kinetic roughen- 
ing are summarized in Figs. [IHa),[T3Ib), andQjJc) which 
shows the variation of global scaling exponents a gi /3, z g 
and local exponents ati, /3* with G. With an increase 
in G from 10 to 10 4 , we found (3 ~ 0.50, a g decreased 
from 0.87(1) to 0.73(1), and z g decreased from 1.71(1) to 
1.38(2). On the local length scale, an increase in G did 
not effect a noticeable change in ol\ (a/ ss 0.46) and /3* 
was observed to decrease from 0.23(4) to 0.18(8). 
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FIG. 14: Variation of exponents a g , ft, z g , ati, and /3» with 
G. The global exponents were calculated after t — 5000 and 
the data were averaged over 1.8 x 10 3 independent runs. The 
local exponents were calculated using L = 2000 and t = 1000. 
The data were averaged over 200 runs and the error in the 
exponents were obtained from the curve- fit. 



IV. CONCLUSIONS 

We have performed 1+1D Monte Carlo simulation of 
VDP growth process by considering free-monomer de- 
position, free-monomer diffusion, polymer initiation, ex- 
tension and polymer merger. The ratio G of the free- 
monomer diffusion coefficient D to the deposition rate 
F was found to have a strong influence on the film's 
growth morphology. The growth rate R(G) of the poly- 
mer film was found to increase monotonically with an 
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increase in G. This is due to the consequence of an up- 
per diffusion flux of free-monomers. The detailed anal- 
ysis of the surface morphology indicated the presence of 
very different scaling behavior at global and local length 
scales. The kinetic roughening study of film interface 
indicates anomalous scaling and multiscaling. With an 
increase in G from 10 to 10 4 , the global growth expo- 
nent (3 « 0.50 was found to be invariant, whereas the 
global roughness exponent a g decreased from 0.87(1) to 
0.73(1) along with a corresponding decrease in the global 
dynamic exponent z g from 1.71(1) to 1.38(2). The global 
scaling exponents were found to follow the dynamic scal- 
ing hypothesis with z g — a g //3 for various G. With an 
increase in G from 10 to 10 4 , the average local roughness 
exponent ai remained close to 0.46 with ai ^ a g , this 
observation is unlike the ones obtained in self-affine sur- 
faces [TH, H|| . The anomalous growth exponent /3* was 
also found to decrease from 0.23(4) to 0.18(8) with an 
increase in G. Even though our model is in 1+1D as 
compared to the 2+1D experiments, our estimates of on 
and /3* are close to the experimental findings of a 0.5 
to 0.7 and [3 = 0.25 ± 0.03 obtained from the AFM stud- 
ies of linear PA-N films grown by VDP [HI, [5l|. The 



similarity between the experimental and simulational es- 
timates appears to be a coincidence since the dimension 
of the two systems are totally different. We also did not 
observe the changes in the dynamic roughening behavior 
reported by Ref. [23jj, perhaps due to the limitations of 
our current simulation model in considering the effect of 
free-monomer diffusion only. This makes us believe that 
the kinetic roughening of the polymer films is sensitive 
to the specific molecular-level interactions, relaxations of 
polymer chains through inter-polymer interactions, and 
the intrinsic nature of polymerization process that need 
to be accounted for in the future simulations. 
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